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ABSTRACT 

Context. The stability properties of toroidal magnetic flux tubes are relevant for the storage and emergence of magnetic fields in the 
convective envelope of cool stars. In addition to buoyancy- and magnetic tension-driven instabilities, flux tubes are also susceptible to 
an instability induced by the hydrodynamic drag force. 

Aims. Following our investigation of the basic instability mechanism in the case of straight flux tubes, we now investigate the stability 
properties of magnetic flux rings. The focus lies on the influence of the specific shape and equilibrium condition on the thresholds of 
the friction-induced instability and on their relevance for emerging magnetic flux in solar-like stars. 

Methods. We substitute the hydrodynamic drag force with Stokes law of friction to investigate the linear stability properties of toroidal 
flux tubes in mechanical equilibrium. Analytical instability criteria are derived for axial symmetric perturbations and for flux rings in 
the equatorial plane by analysing the sequence of principal minors of the coefficient matrices of dispersion polynomials. The general 
case of non-equatorial flux rings is investigated numerically by considering flux tubes in the solar overshoot region. 
Results. The friction-induced instability occurs when an eigenmode reverses its direction of propagation due to advection, typically 
from the retrograde to the prograde direction. This reversal requires a certain relative velocity difference between plasma inside the 
flux tube and the environment. Since for flux tubes in mechanical equilibrium the relative velocity difference is determined by the 
equilibrium condition, the instability criterion depends on the location and field strength of the flux ring. The friction-induced insta- 
bility sets in at lower field strengths than buoyancy- and tension-driven instabilities. Its threshold is independent of the strength of 
friction, but the growth rates depend on the strength of the frictional coupling between flux tube and environment. 
Conclusions. The friction-induced instability lowers the critical magnetic field strength beyond which flux tubes are subject to grow- 
ing perturbations. Since its threshold does not depend explicitely on the friction parameter, this mechanism also applies in case of 
the quadratic velocity dependence of the hydrodynamic drag force. Whereas buoyancy- and tension-driven instabilities depend on the 
magnetic field strength alone, the dependence of hydrodynamic drag on the tube diameter gives rise to an additional dependence of 
growth times on the magnetic flux. 

Key words. Magnetic fields - Magnetohydrodynamics (MHD) - Instabilities - Sun: magnetic fields - Stars: magnetic fields 
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^. Introduction 

X 

^ The magnetic field observed on the surfaces of cool, solar- The stability properties of magnetic flux tubes are rele- 
_C3 like stars is assumed to originate in the bottom of the stel- vant for the amplification, storage, and emergence of mag- 
lar convection zone. It is amplified by the shear flow in the netic fields. An efficient generation of magnetic flux requires 
tachocline and stored in the stably stratified overshoot re- the field to remain within the convection zone for time peri- 
gion at the interface to the radiative core. The toroidal mag- ods which are comparable with or longer than the amplifica- 
netic fields are subject to perturbations through overshoot- tion timescale of the (here, unspecified) dynamo mechanism, 
ing gas plumes penetrating from the convection zone above. Instability mechanisms limit the storage time and lead to their 
Beyond a critical field strength, unstable perturbations lead to leakage of magnetic fields from the overshoot region with field- 
growing magnetic flux loops, which rise through the convec- dependent growth times. The main driving mechanism of un- 
tion zone an d eventually emerge a t the stellar surface (e.g. stable, non-axisymmetric perturbations is magnetic buoyancy^ 
ISpruit & vaiTB allegooijen 1982; Schiissler et alJ [1996). In the which gives rise to undulatory, Parker-type instabilities ( Parked 
solar oversho ot region, critical field streng t hs are of the order |l955l) . If the flux tube dynamics is dominated by the magnetic 
of 10^ G (e.g. lFerriz-Mas & Schussleilll993Lll995h . Simulations tensioi i force, an a xial-symmetric poleward slip instability can 
of rising flux tubes with such initial field strengths yield emer- occur jSpruit & va n Ballegooiien 1982). 

gence properties which are consistent with observed properties The hydrodynamic drag force reduces the relative motion of 

of sunspots such as their eruption latitudes, tilt angles, proper a flux tube perpendicular to its environment, similar to the case 

motions, and asymmetries between the le ading and the foUow- of a solid cylinder immersed in an external flow. The intrinsi- 

ing component of bipolar spot groups jD'Silva & Choudhuril cally dissipative interaction with the environment can cause un- 

der certain conditions the onset of flux tube instabilities (e.g. 

lRvutoval[i988t iJoarder et alJll997 ). A friction-induced instabil- 

Send offprint requests to: V. Holzwarth ity occurs if the relative flow velocity along the magnetic flux 
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tube is higher than a critical value. In paper II of this series 
dHolzwarth et al.ll2007h . we have investigated friction-induced 
instabiUties of straight, horizontal flux tubes with parallel flows 
to elucidate the basic instability mechanism. In that case, the ve- 
locity difference between the internal and external plasma was a 
free parameter. Here we consider magnetic flux rings in a strat- 
ified environment to investigate the influence of the toroidal ge- 
ometry and of the equilibrium condition on the stability prop- 
erties. Both aspects affect the velocity difference between the 
internal and external plasma. 

In Sect. |2] we describe the linearisation of the equations of 
motion in the framework of the thin flux tube-model, including a 
frictional interaction with the environment. In the linear stability 
analysis in Sect. [3] analytical instability criteria are derived for 
some special cases. The general case is analysed numerically 
on the basis of flux rings located in the solar overshoot region. 
Section |4] contains the discussion of the results and Sect. |5] the 
conclusions. 



2. Model description 

2.1. Thin flux tube approximation 

The stability analysis is c arried out in t he framework of the thin 
flux tube approximation ("S pruitlfTQSl I). The radius of the flux 
tube is small compared with the radius of curvature, the wave- 
length of a perturbation, and the pressure scale height. Cross- 
sectional variations are taken to be negligible. All quantities are 
given by their values on the tube axis and are functions of time, f, 
and arc length, s, only. The flux tube is assumed to always retain 
a circular cross section. 

The dynamics of the magnetic flux tube is determined 
through the equation of motion 

p—u = -V\p + —\+pg+n—K+t — — +2p{uxSl)+foAl) 
at \ on I An os ott 

where u is the plasma velocity, B is the magnetic field, p is the 
density, p is the gas pressure, k is the curvature of the tube, and £i 
is the constant rotation vector of the star. For most but the fastest 
rotating stars the contribution of the centrifugal acceleration to 
the effective gravitation is small. We therefore take g ^ -gSr, 
with Cr being the radial unit vector in a spherical coordinate sys- 
tem and g {- const) being the local gravitation. The adopted co- 
ordinate system is the tripod co-moving with the magnetic flux 
tube (Frenet basis), consisting of the tangential vector t, the nor- 
mal vector n, and the binormal vector b. The tangential vector 
defines the axis of the flux tube, implying B - Bt. 



2.2. Drag force 

In the limit of infinite conductivity, plasma cannot cross the mag- 
netic surface of a flux tube. Consequently, a motion of the flux 
tube perpendicular to its tube axis implies a flow of the external 
plasma around the tube. Analogous to a cylinder immersed in a 
non-ideal streaming fluid, this flow gives rise to a hydrodynamic 
drag force 



f - I 

/D - -Pe Mrel,± Mrel,± 

na ' ' 



(2) 



where pe is the density of the external plasma, a is the tube ra- 
dius, Cd is the (empirical) drag coefficient, and Mi-ei.i is the com- 
ponent perpendicular to the tube axis of the velocity difference 
between the internal and the external plasma. In the reference 



frame co-rotating with the star, the external plasma is at rest and 
the relative perpendicular flow velocity thus given by 



Miel,j_ - U - (U ■ t)t 



(3) 



In stationary equilibrium, the plasma flow inside the flux tube is 
parallel to the tube axis. Owing to its quadratic dependence on 
the velocity, the hydrodynamic drag force drops out of the lin - 
earised equation of motion (e.g. Ferriz-Mas & Schiissleij[T993h . 
To investigate the principal influence of friction on the stabil- 
ity properties of magnetic flux tubes, we consider the case of a 
linear velocity dependence according to Stokes's law: 

/d fst = -PeffMi-el.i ■ (4) 

The friction parameter, a > 0, is taken to be constant. 
2.2.1. External stratification 

Disregarding any differential rotation or meridional circulation, 
the stellar stratification in the co-rotating reference frame is de- 
termined through 



(5) 



where and pe are the external gas pressure and density, re- 
spectively. The flux tube is in lateral pressure balance with its 
field-free environment: 



P + — 



(6) 



Subtracting Eq. Q from Eq. ([T]i yields 



d Ap B"- \ d B^ Pe 

—u = — g^iS + riK- — + t- — — + 2{uxil) au^, (7) 

dt p 4np p OS on p 

where Ap = p - pe is the density contrast. Equation ^ is the 
starting point for the linear stability analysis. 



2.3. Linearised equation of motion 
2.3.1. Equilibrium condition 

We consider magnetic flux rings in mechanical equilibrium 
located parallel to the equatorial plane. A detailed descrip- 
tion of this equilibrium and how it is obtained is given by 
Moreno-Insertis et al. ( 1992) and CaUgari et al. (1995). 

In the following, the index '0' indicates equilibrium values. 
Since the equilibrium flux ring is axially symmetric, the tangen- 
tial component of Eq. (|7]l vanishes. From the binormal compo- 
nent follows that the density contrast is nil, because buoyancy 
provides the only force component parallel to the rotation axis. 
With Apo = the normal component of Eq. (|7]) yields 



ul + 2M()Ve,0 - c\fl = , 



(8) 



where ca.o = Bq/ ^J4■npo is the Alfven velocity and Ve,o - Q/kq 
is the rotation velocity of the tube's environment. Equation (O 
determines the velocity. 



Mo - Ve.O 



1 + 



ca,o 



- 1 



> Ve ,0 , 



(9) 



of the plasma inside the flux ring which is required to balance 
the magnetic tension force through Coriolis and inertial forces. 
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This internal plasma flow is in prograde direction. We shall refer 
to the dimensionless quantity 



The coefficient matrices are 



ca,o ca,o fo 



ca,o 



ca,o 



Q7?() Qro cos/l{) 



(10) 



as the equilibrium parameter of the flux tube, where Rq is the 
equilibrium radius of the flux ring in cylindrical coordinates, and 
ro and its radius and latitude, respectively, in spherical coor- 
dinates. 



2.3.2. First-order perturbations 

The linear stability properties are determined through the first- 
order perturbation terms of Eq. (|7]i. In the co-moving Frenet ba- 
sis, the tangential, normal, and binormal components are 

Id \ (1 d , 

1— m| ■'()- -t-^ + KQCj^Qino- ti) -2D.(no-Ui) =(1Q) 



[dt 



'1 



d \ An 
I —uj - no- — 



[dt 



p 55 Stt/j 

gOgn +2n(to-Ui) + a (Hrel,±,l ■ «o) 



d \ Ap 
— M I bo 

P 



\dt 



■KoCj,J-+2—-—\ ^1(2) 

\Kq «o Pi) I 



gOgb + KoC^^f (no • bi) H- Q;(Hrel,±,l ' *o) =(KB) 



with g„ = (g ■ no) /go and g,, ^ (g ■ bo) /go- For magnetic flux 
rings perpendicular to the stellar rotation axis it is to - e^, no - 
-en, and bo = e-. The first-order perturbations of all quantities 
(indicated by the index '1') depend linearly on the Lagrangian 
displacement vector ^ - {^,,^„,^hf , and are summarised in 
Append ix [A] Their detailed derivation can be f ound, for exam - 
ple, in Perriz-Mas & Schiissled(ll993Lll995h and lSchmittl(ll998h . 

Except for the drag force, which comes from the perturbation 
of the exterior, the velocity of the external plasma is taken to be 
unchanged by the perturbation of the flux tube, that is, Me,i = 0. 
The perturbations of the internal flow velocity and of the tube's 
axis given by Eqs. ( |A.5t and dA.lt , respectively, yield the per- 
turbation of the relative perpendicular flow velocity: 



«rel,±,l - Ml - ("l ■ to) to - Uot\ - ^njUQ + ^b,tbo 



(14) 



Derivatives with respect to time and arc length are indicated by 
the indices f and s, respectively. At the bottom of stellar con- 
vection zones, the plasma pressure is typically much higher than 
the magnetic pressure. The stability analysis is therefore carried 
out in the approximation y6 = %npolB^^ » 1. In the limit of 
an infinite radius of curvature and vanishing stellar rotation (i.e. 
/<• — > 0, Q — > 0), Eqs. (fTTTi-dmi reduce to the case of a horizontal 
flux tube in a stratified atmosphere, which was investigated in 
iPaper 11 Note that the present investigation dis regards the pos- 
sible influence of enhanced inertia (factor jj. in Paper II), owing 
to i ts conceptual difficulties in t he case of curved thin flux tubes 
(cf. lMoreno-Insertis et al.|[T996l and references therein). 

2.3.3. Non-dimensional form 

The linearised equations of motion are transferred to non- 
dimensional form by multiplying Eqs. (fTni-(fT3Tl with the 
timescale t = 1 which yields 

T't„+Ms,-^,s,+M,s^ts.+M,T^,,+Ms-ls+M^^ = .(15) 



Ms, = 
Ms, = 
M, = 

Ms = 
where 



2Ro£, 
-2Ro6 

2 (Ro + l) + aX 
-2/?oJ/+-^-(g„J/ + gbZ) 

r/o 



,~,2 



(0 \ 

gngv^ gv^gh 
gng 



rf_ 

yfo 



(0 \ 

2gn gh 

gi, 



(16) 
(17) 
(18) 

(19) 
(20) 
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lo 





oj 




[ 1 





j 



(21) 



£ is the unity matrix. 



Ro 



uoKo r. , 



is the Rossby number, 
fo - Hp Ko 



(22) 



(23) 



is the ratio between the local pressure scale height and the radius 
of curvature, and a - Tais the dimensionless friction parameter 
The (dimensionless) magnetic Brunt-Vaisala frequency. 



,~,2 



2, ,2 



= T W- 



2/o' 



/36 



2 






7 




2 



(24) 



is a measure for axisymmetric buoyancy oscillations of mag- 
netic flux rings. Owing to the stabilising effect of the magnetic 
field, its value is higher than the non-magnetic Brunt-Vaisala 
frequency. 



(25) 



Ko 



Ko 



which depends on the superadiabaticity 6 -V - Vad and quanti- 
fies the stability of the stratification against convective turnover. 
Positive (negative) values of w^gy and A^^ indicate stable (un- 
stable) perturbations of the plasma in the presence and absence 
of a magnetic field, respectively. 

2.3.4. Solar reference case 

For the numerical evaluation of linear stability properties, we 
shall refer to a magnetic flux ring located in the middle of the so- 
lar overshoot region. At radius ro = 5.07 ■ 10"^ cm, the local den- 
sity is po = 0.15 g/cm-', the gravitation is go - 5.06 ■ 10"* cm/s^, 
the pressure scale heigh is Hp = 5.52 ■ 10'* cm, and the su- 
peradiabaticity is ^0 = -10"*. The solar rotation rate is taken 
to be = 2.69 ■ 10''' s"' (i.e. P,ot = 27 d) and the ratio of 
specific heats y = 5/3. The local Brunt-Vaisala frequency is 
A^o = 3 ■ 10"^-' andM) = N/Q. = 1.11. 
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3. Linear stability analysis 

3.1. Method based on sequence of determinants 

We consider perturbations in the form of discrete harmonic func- 
tions. The exponential ansatz 

^ - i exp [/ (niKs - Li)Q. t)] (26) 
transforms Eq. ( fTsT i into the homogeneous algebraic system 

[w^ -majMs, + Mss + HhMi - iniMs - M.^\ 1 = 0, (27) 

where oj - + ioji is the (dimensionless) complex eigenfre- 
quency and m is the azimuthal wave number; owing to the 2jt- 
periodicity of the flux ring, the latter is an integral number. The 
determinant of the coefficient matrix of Eq. dZTl i defines a 6th- 
degree dispersion polynomial in tl>, which determines the stabil- 
ity properties of the magnetic flux ring: eigenmodes with d); < 
are (asymptotically) stable, while those with dl>i > are unsta- 
ble. The coefficients of the dispersion polynomial are given in 
Appendix lB.il 

An explicit calculation of eigenfrequencies is not essential 
for the determination of the stability criteria of a dynamical sys- 
tem, which can be inferred from the coefficients of the disper- 
sion polynomial by using theorems from stability and control 
theory. In our case, an instability implies a root of the disper- 
sion equation in the upper half of the complex plane. We there- 
fore use a ge neralised version of the Routh-Hurwitz theorem 
( lMardenlll966t , Sect. 39): Assume the monic complex polyno- 
mial fiz) = C0 + C1Z+ ■ ■ ■ + c„_iz""' +z" i^O for z real, and define 
the quantities 



with the real coefficients C2 and co given in Eqs. dB.Sb and ( IB.9b . 
respectively. The elements of the 'V-sequence are: 



A7. = 



1 I 
-Di = - 

2 2 



— £>2 

4 ^ 



1 2ia 
1 —2ia 



2a 



2ia C2 ■ 
-2 ia C2 ■ 



2ia 



iac2 
-i a C2 

C2 



-2 ia C2 - 



-2a^ [c2-2a^) 



A, = 



1 


2 i a 


C2 - 


iac2 


Co 





1 


—2 ia 


~2 

C2 — a 


—i a C2 


CO 








1 


2 ia 


~2 

C2 — a 


iac2 


CO 





1 


—2 ia 


~2 

C2 - a 


-iac2 


Co 








1 


2ia 


~2 

C2 - a 


iac2 








1 


—2ia 


~2 

C2 — a 


-iac2 



= -a^ (c2 - 2 a^) (c\ - 2 5-^ C2 - 4 cq) 



A4 = 



Da 



16 

a'^ Co {c\ ■ 



16 



2 a C2 - 4 Co j 



(31) 



(32) 



(33) 



(34) 



According to A2, A3, and A4, instabilities occur if C2 > 2a^ or 
Cj - 4co - 25'^ C2 < or Co < 0, respectively. 
The A2-criterion implies the condition 



through the determinants Dk of the 2k x 2A:-principal minors (i.e. 
the first 2k elements in the first 2k rows and columns) of the 
square 2n x 2«-matri)|3 
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C„-l 


C„-2 




Co 








. . . '\ 
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C„-l 


C„-2 




Co 
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C„-l 


C„-2 




Co 
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C„-l 


C„-2 




Co 












1 


C„-l 


C„-2 




Co 










1 


C„-l 


C„-2 




Co 
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C„-I 


C„-2 




Co 


v ■ ■ ■ 








1 


C„-I 


C„-2 




Co . 



(29) 



Then if A^ for k - 1, 2, . . . , n, the number of zeros of /(z) in 
the upper half of the complex plane is equal to the number sign 
changes in the sequence 'V = {I, Ai, A2, . . . , A„). 

The method is used in the following sections, which focus 
on special cases of flux tube perturbations. 

3.2. Axisymmetric perturbations 

The flux tube equilibrium is indifferent with respect to axially 
symmetric perturbations in the azimuthal direction. Such pertur- 
bations are marginally stable and the associated eigenfrequen- 
cies zero. The remaining eigenfrequencies are determined by the 
quartic dispersion polynomial 



+ 2 i a 0/ + (c2 — a^) ci/ + i a C2 ci) + cq — , 



(30) 



A bar over a complex quantity denotes its complex-conjugate value. 



(28) '^MBV<^^^-2j/^.(/^o.2)----^a 



^2 



(35) 



In the frictionless case, or = 0, th is criterion describes a 
monotonic buoyancy-driven instability dFerriz-Mas & Schiissleii 
Il995 ). Here, friction has a stabilising effect and reduces the crit- 
ical magnetic Brunt- Vaisala frequency below which the instabil- 
ity sets in. 

In our analysis of th e A3 -criterion, we follow the approach 
of iFerriz-Mas & Schiis sler (1995) and combine the coefficients 
Co and C2 to eliminate tt'^gy- resulting expression is used 
to eliminate co from the A3 -criterion, which yields a quadratic 
expression in C2. If the discriminant. 



8l 



^ ^ a'' - 16^ (Ro +!)'&■ 

4 n2 



2 ~2 



8l 



—Ro (Ro + 2) -4— (Ro + if 

yfo g± 



(36) 



of this quadratic expression is positive and the magnetic Brunt- 
Vaisala frequency in the range 



~ 2 ~ 2 
a)^_ < cl). 



~ 2 

MBV ^ '^c,+ 



where 



~2 



'4 



2gn 

r/o 



(37) 



(38) 



then A3 is negative. In the frictionless case this criterion de- 
scribes a buoyancy-driven overstabili ty which occurs if the 
enviro nment is rotating differentially dFerriz-Mas & Schiissleii 
Il995h . 
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The A4-criterion applies to flux tubes outside the equatorial 
plane only and implies 



4\r/o/ i+if ' 



(39) 



This instability is driven by the tension of magnetic field lines 
and causes a monotonia movement of the flu x ring to higher lat- 
itudes (cf. ISpruit & van Ballegooiienlll982h . In contrast to the 
A2- and A3 -criteria, the threshold of this 'poleward slip instabil- 
ity' depends neither on the friction parameter nor on the local 
pressure scale height; the latter is indicative for buoyancy-driven 
instabilities. 

The poleward sUp instability is the governing mechanism re- 
garding axially symmetric perturbations of flux rings outside the 
equatorial plane. This is shown by verifying that the sign change 
in the 'V-sequence first occurs between its last two elements. At 
the A4-threshold, defined by equality of both sides of relation 
( [39] l. that is cq = 0, we have 



1 

r/o 
1 

r/o 



Tf 



4 (Ro + 1) 

4 {Ro+ \ f ' \r/o 



2 



+ 2\ — -2\if - A 



< 



4(1-1- if) 
and, by applying these relations, 

[cl - 4co - 2a\2\^^^ = [cl - 2a\2\^^^ 

= [c2]co=o [c2 - 2a\^^^ > . 



(40) 



(41) 



A magnetic flux ring which becomes subject to the poleward slip 
instability is therefore still insusceptible to buoyancy-driven in- 
stability mechanisms. This applies to flux rings outside the equa- 
torial plane only, where the absolute term of the dispersion equa- 
tion is not zero per se. 

In the equatorial plane, buoyancy-driven instabilities occur 
if C2 > 0. Since the discriminant dc^ = is positive, the A3- 
criterion covers the range 



2(llL_2U2-4-2a2<£;2 

\r/o 



MBv<2(^-2l,^-4. 



(42) 



and therefore extends the instabiUty criterion ( |35] | to higher mag- 
netic Brunt- Vaisala frequencie s. In analogy with the results of 
iFerriz-Mas & Schiisslej (Il995h for the frictionless case, we ex- 
pect that in the presence of diff'erential rotation the Aa-criterion 
describes a buoyancy-driven overstability. In the present case of 
solid body rotation, explicit calculations of the eigenfrequencies 
show that (within computational accuracy) axially symmetric in- 
stabilities are monotonia. 

In the determination of the critical equilibrium parameter, 
beyond which magnetic flux rings are subject to unstable axi- 
ally symmetric perturbations, the dependence of the magnetic 
Brunt-Vaisala frequency on the field strength must be taken into 
account. To this end, we use Eq. jTM to substitute w^gy in Eqs. 
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10^ 
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Fig. 1. Thresholds of the poleward slip instability and axially 
symmetric equatorial instability for different length scale ratios 
/o = Hp kq. Flux tube equilibria with parameters located above 
(below) a curve are insusceptible (susceptible) to the respective 
instabihty. 



( [39] l and ( |42] |. The equation for the threshold of the poleward slip 
instability is 



(3-2.)^-|„V 



+ (4- ly) 



1 



7,2 +4^r2 ^0.(43) 



This instability occurs in a stably stratified environment (i.e. 
N~ > 0) if the equilibrium parameter is sufficiently large (Fig. 
[TJ. In contrast, a flux ring located in a convectively unstable en- 
vironment is insusceptible to this instability mechanism if the 
local Brunt-Vaisala frequency is higher than 



mm 



3-r + i 



-1.6- 10-2 /./o 



(44) 



and the equilibrium parameter is within a range of values delim- 
ited by the solutions of Eq. ( |43] ); an upper boundary of this range 
is 



j](N^ = 0) 



4-2r 



V2r-3 



V2. 



(45) 



The threshold of the buoyancy-driven instability in the equatorial 
plane is determined by 



(46) 



In a stably stratified environment a flux ring is susceptible to this 
instability mechanism provided that the equilibrium parameter 
is sufficiently large and the ratio between pressure scale height 
and radius of curvature is in the range 



1- V2r^ ^ 1+ V2r^ 

< ./() < ^n- 



4y 



4y 



(47) 



If the length scale ratio is outside this range, the axisymmetric 
equatorial instability can only occur in a convectively unstable 
environment (Fig.[T] thin lines). With 7 = 5/3 and = 1, the 
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range is 0.063 < /o < 0.24. In the overshoot region of the Sun it 
is (at the equator) /o a; 0.1. 

Both poleward slip and equatorial instabilities are not caused 
by friction in the first place and occur in the frictionless case 
as well. There are no specific friction-induced instabilities re- 
garding axially symmetric perturbations of magnetic flux rings, 
though friction may affect the occurrence of overstabilities in the 
presence of differential rotation. The absence of friction -induced 
instabilities is consistent with our finding in iPaper III that their 
occurrence depends on the relation between flow and phase ve- 
locity. Since axially symmetric perturbations do not yield propa- 
gating wave signals the method of principal minors recovers the 
(modified) instability mechanisms which also occur in the fric- 
tionless case. 



3.3. Flux tubes in the equatorial plane 

We now consider flux rings in the equatorial plane and derive an- 
alytical criteria for friction-induced instabilities with azimuthal 
wave numbers m > 1. In the equatorial plane binormal pertur- 
bations are decoupled from perturbations within the osculating 
plane, so that the dispersion equation factorises into a quadratic 
and a quartic polynomial. 

The eigenfrequencies of binormal perturbations, 



_4(^)\=K-3).2(i.)',. 

In Eqs. (I50l)-(l53l). the friction parameter occurs as a multiplica- 
tive factor only: oc a''. The instability thresholds therefore 
depend on the presence of friction only, but not on its strength. 

It is sufficient to analyse the lowest-order eigenmode, be- 
cause A2, A3, and A4 increase with the azimuthal wave number: 



dm 



> 



Ak(m\) < Ak(m2) for < nii < 1112 (56) 



The increase of A3 with m is not obvious, since it depends on F\ . 
Regarding F\ as a quadratic function in for example, its roots 
ga,± are complex-conjugate, so that for real-valued ga the sign of 
F\ is unique. Since the function is positive for, say, g„ - 0, it 
foflows Fx > 0. 

An instability is indicated by A2 < 0, that is. 



2 (Ro + 2) 



r/o 



2\Ro-\ 



or by A3 < 0, implying 
F2-2R0 Fi 



r2 



< 



(57) 



(58) 



■1 



m Ro — i — ± 
2 



i a m Ro 

4 



(48) 



reduce in the frictionless case to m (Ro + rf). The corresponding 
phase velocities Vph = mo ± ca imply perturbations which propa- 
gate with the Alfven velocity in the prograde and the retrograde 
direction relative to the internal plasma flow. In the frictional 
case, a binormal eigenmode is unstable if 



3 



if- — i a m Ro — — 



a 
> — , 
2 



(49) 



which requires flow velocities mq > ca. This requirement is 
never fulfilled, because the mechanical equilibrium permits sub- 
Alfvenic flow velocities only. Binormal perturbations of flux 
rings in the equatorial plane are thus stable both in the presence 
and in the absence of friction. 

The eigenfrequencies of perturbations in the osculating plane 
are determined by a monic quartic polynomial, whose coeffi- 
cients are given by Eqs. ( |B.10l )-( |B.13t in Appendix IB. 3 1 The 
associated 'V-sequence is 



Ai 

A2 



«2 



A[\+if-) + 2m'Ro-2-^ri'' +gl: 
A3 = 1]^ [(2 Ro + gl ^^3^) Fi - F2] 



MBV 



(50) 
(51) 
(52) 



A4 = a'^m^i]'* 



6n 1 2 /I fin 2 a 

— r] -4 — T] -4 

r/o / r/o 



4 (m2 - 1) Ro^ - j i^U2gl w^Bv 



(53) 



with 



F, = 4 (47,2 - 2Ro + 1) + 4^ [Ro - 2if-) + j 77^ (54) 

F2 = l6il+Tf){2Ro-\)-%—RoiATf- + 2,Ro + %) 
^ ' r/o ^ ' 



The A4-instability occurs for 

v2 



1 / 1 



(59) 



We shall show that the first sign change in the sequence 'V occurs 
at its end, so that at the threshold of the A4-instability both A2 
and A3 are still positive; from a > follows Ai > 0. After 
substitution of ^I'^bv ^1- 'EB with the right side of Eq. ( |59l ) 
we get 



52 



[A2]c„=o ^^iRo + 2) 



77^+4 2-^ /?o + 4 

/or/ I r/o' 



.(60) 



Following the same argumentation as after Eq. ( 1561 ) for function 
Fx, it can be shown that the term in square brackets is positive. 
Therefore, A2 > 0. For A3 we find 



{A^\,^^Q^—Tf{Ro + 2) 



r/o/ ' Ir/o' ' 



(61) 



which is non-negative. Relation ( |59] | is therefore the instability 
criterion for non-axial symmetric perturbations of magnetic flux 
rings in the equatorial plane. 

Taking the field dependence of the magnetic Brunt-Vaisala 
frequency into account, the threshold of the A4-instability is 
given by 



1/ 1 



= 2\^J Ro{Ro + 2){Ro + y) . 



(62) 



Its dependence on the equilibrium parameter is shown in Fig. 
|2] together with the thresholds of the A2- and A3-criteria. In the 
equatorial plane, magnetic flux rings are unstable against non- 
axial symmetric perturbations if they are located in a convec- 
tively unstable stratification. If the environment is stably strati- 
fied, the instability occurs provided the equilibrium parameter is 
sufficiently large. 

Since there is only one sign change in the 'V-sequence, only 
one eigenmode is unstable yielding growing perturbations. A 
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Fig. 2. Thresholds of the friction-induced instability (solid line) 
and of buoyancy-driven instabilities (broken lines) for non-axial 
symmetric perturbations (for /o - 0.1). The dotted, horizontal 
line indicates the stratification parameter .^V^ = 1 .23 in the mid- 
dle of the solar overshoot region. 




0.2 0.3 0.4 0.5 



Fig. 3. Growth rates, of magnetic flux tubes in the 

equatorial plane. The threshold to the left is determined by the 
friction-induced instability with azimuthal wave number m - 1 . 
In the region enclosed by dotted lines, m - 2 modes have 
the highest growth rates. The dimensionless Brunt-Vaisala fre- 
quency is N^/Q? = 1.23 (cf. Figllli and/o = 0.109. 

comparison with Eq. (l24l i shows that Eq. ( |62| | is independent 
of the pressure scale height, so that the underlying instability is 
not caused by buoyancy. It does not exist in the frictionless case 
and is therefore specifically induced by the frictional interaction 
between the magnetic flux tube and its environment. The driving 
mechanis m of the fr iction-induced instability has been described 
in detail in 'Paper 11*. 

The basic condition which discriminates between stable and 
unstable flux tube equilibria is the A4-criterion with azimuthal 
wave number m = 1 . Beyond that threshold, higher-order eigen- 
modes can be unstable with growth rates higher than those of 
the m - 1 mode. The determination of the fastest growing per- 
turbations requires the explicit calculation of all eigenfrequen- 
cies. Keeping the equilibrium parameter constant, growth rates 
decrease if friction is strong (Fig. [3]). This is caused by the dissi- 



pative nature of friction, which opposes relative plasma motions 
perpendicular to the tube axis. If friction is weak, the dissipative 
effect is outbalanced by the energy transfer from the relative mo- 
tion of internal and external plasma through frictional coupling. 
Further to friction, the growth times of unstable perturbations 
still depend on magnetic buoyancy and tension. Growing pertur- 
bations of flux tubes which in the frictionless case are subject 
to Parker-type instabilities experience in the frictional case still 
a strong magnetic buoyancy which, gives rise to relative high 
growth rates (cf. Fig. [3] lower-right corner). 

An important aspect of our approach using the principal mi- 
nors of the coefficient matrix D is that, for a polynomial of de- 
gree n with 3 (cq) = 0, the last in the sequence depends lin- 
early on the absolute term of the dispersion polynomial: A„ oc cq. 
This can be seen by comparing Eq. ( l53T l with Eq. ( IB. 131 ) and also 
in Eq. ( [34] l. A sign change of the absolute term of the dispersion 
polynomial thus indicates a friction-induced instability with, at 
least, one eigenfrequency being zero at the th reshold cq = 0. 
This result agrees with our finding in iPaper iS that the onset of 
a friction-induced instability is characterised by the reversal of 
the propagation direction of a wave mode from retrograde, with 
(w) < 0, to prograde, with % (fij) > 0. 

3.4. General case 

In case of magnetic flux rings parallel to but outside the equato- 
rial plane, gravitation causes a coupling of perturbations parallel 
and perpendicular to the rotation axis. The eigenfrequencies are 
determined by a 6th-degree dispersion polynomial, whose co- 
efficients are given in Appendix lB.il The method based on the 
principal minors is still applicable, but the analytical expressions 
for Ak are cumbersome and not of much practical use. Guided by 
our previous results, we take the reversal of the propagation di- 
rection of an eigenmode as the criterion for the friction-induced 
instability, implying A(, oc -co < 0. The resulting instability 
criterion is formally identical with condition ( |59] l. 

The latitude-dependence of the instability threshold, Eq. 
( |62] l, is mediated through the dependence of Rossby number 
and equilibrium parameter on the curvature of the equilibrium 
flux tube. Consider, for example, similar magnetic flux rings at 
a given depth of a spherically symmetric star. Since the Brunt- 
Vaisala frequency is the same for each flux tube at different lat- 
itudes, the critical flow velocity decreases with increasing lat- 
itude: Mocrit cos /!(). Owing to the mechanical equilibrium 
condition (|9]l, this implies that at higher latitudes the friction- 
induced instability sets in at lower field strengths (Fig.|4|i. If the 
latitude dependence of all parameters is taken into account, the 
threshold is determined through the bi-cubic polynomial 

T/f - (r - 2) y cos^ A7]'^^+4-(y-l)y-Bu cos^ ^t?^ 

- cos^ ^ = . (63) 

The index V of the equilibrium parameter rjr = CA.o/Oro indi- 
cates that this quantity is a function of the equilibrium depth, 
ro, since all latitude dependencies are now explicitely taken 
into account through the cosine-terms in the coefficients. For 
/i = 0, Eq. ( l63b is equivalent to Eq. ( |62] |. The Burger num- 
ber Bu = (NHp/Q.rQ)^ is a measure for the ratio of buoyancy 
and Coriolis force, with low (high) values indicating dynamics 
dominated by rotation (stratification). Fig. |4] shows a compar- 
ison of the instability threshold for different Burger numbers, 
corresponding to, for example, different stellar rotation rates or 
Brunt-Vaisala frequencies. 

- The minus sign is caused by the factor (;/2)* in Eq. l l28t . 
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Fig. 5. Growth rates of Parker-type and friction-induced instabilities of magnetic flux rings in the middle of the solar overshoot 
region. Shaded areas indicate unstable equilibria, with dotted lines marking boundaries of Parker-type instabilities in the frictionless 
case (panel a). Growth rates first increase with friction (panel b), then reach a maximum (panel c), and decrease for very strong 
friction (panel d). Arrows in panel b indicate the latitudes shown in Fig.|7] 
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Fig. 4. Latitude dependence of the friction-induced instability 
threshold for different Burger numbers; in the solar reference 
case, it is Bu - (NH^IQ.rQf = 1.47 ■ 10"^. Note that the equato- 
rial value 77r_eq incrcascs with the Brunt-Vaisala frequency of the 
stratification (i.e. for high Bu). 



We numerically investigate the friction-induced instability 
properties of flux rings outside the equatorial plane on the basis 
of the solar reference case (Fig. |5]l. The calculation of the 'V- 
sequence confirms that the A6 -criterion determines the onset of 
friction-induced instabilities (Fig.|6]l. For field strengths beyond 
this threshold, the sign change in the 'V-sequence can occur at 
A5 or A4, but in the parameter range considered here there is 
always only one sign change, indicating a single unstable eigen- 
mode. Given that the friction-induced instability sets in once an 
eigenmode reverses its propagation direction, the retrograde and 
prograde slow wave modes are liable to become unstable since 
the moduli of their phase velocities are smallest. Owing to the 
prograde plasma flow inside the flux tube, the retrograde wave 
is advected against its propagation direction and (the modulus 
of) its phase velocity closest to zero. In most cases it is the ret- 
rograde eigenmode which becomes frictionally unstable once its 
propagation reverses into prograde direction. An exception to 
this rule occurs if the threshold of the friction-induced instability 
coincides with strong driving caused by magnetic buoyancy. In 
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Fig. 6. Sign changes of A4, A5, and As in the solar refer- 
ence case, coiTesponding to the instability diagrams in Fig. |5] 
Tickmarks along each line indicate the downhill direction of the 
A^ -functions. Within the shown range of field strengths and lati- 
tudes, Ai_3 > 0. 



the frictionless case, Parker-type instabilities are characterised 
by a merging of two eigenmodes: the phase velocities of both 
modes are identical and the growth rates of opposite sign (cf. Fig. 
|7] panel a). If the merging coincides with the threshold of the 
friction-induced instability, the phase velocities of both eigen- 
modes are zero (panel c). In that situation, the prograde eigen- 
mode can reverse its propagation to the retrograde direction and 
become frictionally unstable (panel d). 

In contrast to the frictionless case, where magnetic flux 
rings at high latitudes may be stable up to high field strengths, 
flux rings subject to friction are generally unstable if their field 
strength is beyond the threshold given by Eq. ( |62] i. The growth 
rates depend on the dominating driving mechanisms, which in 
the parameter range considered here are frictional coupling and 
magnetic buoyancy; tension-driven poleward slip instabilities 
typically occur at higher field strengths. If friction is weak, 
growth rates are significant only if the flux tube equilibrium is in 
a parameter domain of Parker-type instabilities (Fig.|5] panel b). 
If friction is strong, a differentiation between Parker-stable and 
Parker-unstable domains becomes irrelevant, since the growth 
times are comparable. In case of very strong friction, growth 
rates decrease as flux tube motions are more and more harnessed 
by the environment. 
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Fig. 7. Oscillation frequencies of slow eigenmodes near to the 
threshold of friction-induced instabilities in the solar reference 
case. Each panel corresponds to a latitudinal cut in the sta- 
bility diagram Fig. |5] panel b; the friction parameter is o; = 
10"' s '. Thick lines indicate an unstable eigenmode. In the 
gray shaded regions the instability is predominantly driven by 
magnetic buoyancy, outside by frictional coupling. Since the az- 
imuthal wave number is m = 1, we have Re(a>/Q.) = Vph/vrot- 



4. Discussion 

The criteria for friction-induced instability of straight and 
toroidal magnetic flux tubes are identical. From the equilibrium 
condition, Eq. (|9]l follows 



Ro = 



I- Ml ' 



(64) 



where Ma = mq/ca is the Alfvenic Mach number, which is lim- 
ited in the toroidal case to values Ma < 1. Substitution of the 
Rossby number in the instabiUty criterion, Eq. ( |59b , yields 



1 



ca 



[yHpj 1-Ml 



(65) 



which is the instability criterion of straight, hori zontal flu x tubes 
in the limit of long wavelengths (cf. Eq. (32) in lPaper flh . In the 
limit of an infinite radius of curvature Eq. ( |65] ) becomes 



/3d > 



1 

r 



(66) 



since mq — > for /cq — > 0. Relation ( 166b is the instability cri- 
terion of a straight flux tube in a stratified environment with- 
out a velocity difference betwee n internal and external plasma 
jSpruit & van Ballegooiienlll982l) . 

In contrast to the case of straight flux tubes, in which the 
flow velocity is a free parameter, the relative flow velocity inside 
toroidal flux tubes is fixed by the equilibrium condition. This 
implies an explicit dependence of the friction-induced instabil- 
ity on the curvature and field strength of the flux ring. The in- 
stability criterion, which relates the flow velocity with the phase 
velocity of a perturbation, thus determines a critical magnetic 
field strength (or, more general, Alfven velocity) beyond which 
the friction-induced instability sets in. This critical field strength 
is lower than the threshold of the buoyancy-driven instability. A 
preliminary numerical parameter study indicates that this prop- 
erty applies to a broad range of stellar rotation rates and equi- 
librium depths in the overshoot region. From this, we conjec- 
ture that it is the friction-induced instability which determines 



whether flux ring equilibria in the overshoot region are stable or 
not. 

Virtually all flux tube equilibria beyond the instability 
threshold are unstable. This is important, for example, in the 
case of fast stellar rotation. S tability analyses of flux rings 
in rapidly rotating stars (e .g. iFerriz-Mas & Schiissleii 119951: 
iHoIzwarth & S chiissler 2003) show that, in the frictionless case, 
for high magnetic field strengths isolated regions of stable equi- 
libria exist, in which flux tubes are not subject to the Parker-type 
instability (see Fig. [8] panel a). These field strengths are well be- 
yond the threshold of the friction-induced instability. Therefore, 
if friction is taken into account, these equilibria are unstable with 
high growth rates (Fig.|8] panels c and d). 

The possibility to store magnetic flux in the overshoot region 
depends on the growth rate of the instabilities. An evaluation of 
the growth rates of friction-induced instabilities requires an es- 
timate for the friction parameter a. A comparison of drag force, 
Eq. (I2]i, and friction force, Eq. (|4]i, yields 



Cd 



Mrel,± 



(67) 



The friction parameter of a magnetic flux tube which gives rise to 
a small sunsp ot is estimated thro ugh non-linear numerical sim- 
ulations (see Calig ari et alJll995l for a description of the code). 
With Bo = 10^G,ao = 1000 km (i.e., O ^ 3 ■ lO^' Mx) and 
Cd = 1 we obtain a ~ 10"' s"'. The results of the stability 
analysis in Fig. |5] show that for this value of the friction param- 
eter growth rates are several orders of magnitude lower than the 
rotation rate. Sunspot-producing flux tubes are thus hardly sus- 
ceptible to friction-induced instabilities. 

Thin magnetic flux tubes are weeded out of the overshoot 
region on shorter time scales than thick flux tubes, and should 
therefore occur more often. This agrees with the trend indi- 
cated by umbral areas observed on the Sun, which follow a log- 
normal distribution with small s urface features being more nu- 
merous than large features (e.g. iBaumann & Solankill2005l) . If 
the higher number of small surface features is a consequence of 
higher friction-induced emergence rates, the dynamics and sta- 
bility properties of thin magnetic flux tubes are dominated by 
friction. In that case, the latitudinal distribution of surface fea- 
tures may also depend on the size of the magnetic feature, since 
thinner flux tubes experience a stronger deflection to higher lati- 
tudes by meridional circulation (e.g. Holzwarth et al. 2006). The 
influence of meridional circulation would possibly be diluted by 
the action of convective motions on rising flux tubes. But since 
the former effect is systematic and the latter random, a discern- 
able signature may persist in the latitudinal distribution. Even 
if thin flux tubes do not make it to the surface due to convec- 
tive interaction, their susceptibility to friction-induced instabili- 
ties implies a decrease of magnetic flux in the overshoot region, 
which affects the budget and amplification of magnetic fields in 
the tachocline region. 

The numerical results shown in Sect. 13.41 are based on mag- 
netic flux rings located in the middle of the solar overshoot layer 
Owing to the strong dependence of stability properties on the su- 
peradiabaticity, flux tubes located deeper in the overshoot region 
require significantly higher field strength to become buoyancy 
unstable. The friction-induced instability eases the need for very 
high field strengths, since its threshold is below the critical field 
strength of the buoyancy-driven instability. For example, a flux 
ring located in the equatorial plane close the radiative core, in an 
environment with 5 - -lO"'* (A^ = 3-10"^ s"'), becomes friction- 
ally unstable for field strengths above 3.5 ■ 10^ G, which is well 
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Fig. 8. Like Fig.|5] but for a rapidly rotating star with a rotation period of P^ot = 2d. The Burger number is Bu = 8.04 ■ 10"^. In 
accordance with the results shown in Fig.H] the latitude dependence of the threshold is weaker than in the solar reference case and 
only effective at very high latitudes. 



below the critical field strength of the Parker-type instability at 
6.2 - lO^G. 

The friction-induced instability of toroidal flux tubes is also 
relevant in other astrophysical contexts, such as in circumstellar 
accretion discs. Accretion of mass implies the transport of an- 
gular momentum to larger radii, which requires an efficient cou- 
pling of adjacent disc annuli. The coupling of Keplerian shear 
flow through weak magnetic fields gives rise to the magneto- 
rotational instability (MRI), which generates magnetic turbu- 
lence and efficiently increa ses the viscous couplin g and angu- 
lar momentum transport ( Balbus & HawlevI Il99lb . The MRI 
requires a poloidal magnetic field and is independent of the 
toroidal field component. This axi-symmetric shear instabil- 
ity occurs for vertical wavelengths larger than a critical, field- 
dependent value, though growth rates and threshold are inde- 
pendent of field strength. It is an interesting possibility that 
the friction-induced instability of toroidal fields complements 
the MRI by supporting the radial magnetic coupling of a disc. 
Both instability mechanisms produce their highest growth rates 
on large length scales, yet the friction-induced instability does 
intrinsically not saturat^ at high field strengths as the MRI. 
However, it is unclear what mechanisms aside of Keplerian shear 
can produce the required velocity difference between the plasma 
inside and outside the flux tube. 

ISchramkowskil ( Il996h investigated the dynamics and evolu- 
tion of small flux rings in an accretion disc, including the drag 
force caused by shear flows. If a flux ring is deformed by a strong 
shear flow, variations of the tension force and density contrast 
occur along the flux tube, which produce two loops rising ver- 
tically through the disc until they emerge into the disc corona. 
These investigations confirm the strong influence of the drag 
force on the dynamics of flux tubes with small minor tube radii, 
but they do not consider friction-induced instabilities. A stabil- 
ity analysis including the drag force would be possible, because 
there is a flow component perpendicular to the axis of the equi- 
librium flux tube. Since differential rotation has a significant in- 
fluence on the stability prope rties of magnetic flux rings (e.g. 
iFerriz-Mas & SchussleJ 11991) . it will be necessary to include 
the effect of shear flows and more general equilibrium condi- 



^ Since the critical wave length, beyond which the MRI occurs, in- 
creases with the field strength, it can exceed the vertical extend of the 
disc if the field strength is sufliiciently high. 



tions in the stability analysis of the friction-induced instability 
in the framework of accretion discs. 

It is not clear to which extend the results based on a linear 
velocity dependence of friction carry over to the non-linear case. 
Given the or-independence of the instability thresholds, the ac- 
tual functional form of friction, that is, linear or quadratic, seems 
to be of secondary importance, provided that it is anti-parallel 
to the flow velocity perpendicular to the tube axis. Numerical 
simulations based on the drag force confirm the properties of 
friction-induced instability derived here. A detailed analysis of 
the non-linear case and parameter study is the topic of a forth- 
coming paper in this series. 

5. Conclusion 

The friction-induced instability determines the stability proper- 
ties of flux rings at the bottom of the convection zone, as its 
threshold is at lower field strengths than those of buoyancy- 
and tension-driven instabilities. The stability properties depend 
on the magnetic field strength and on the magnetic flux. The 
threshold also depends on the location of the flux ring, since 
the toroidal geometry causes a connection between the instabil- 
ity criterion and the equilibrium condition. The analytical re- 
sults confirm that a magnetic flux tube subject to friction be- 
comes unstable when a backward wave occurs, that is, when an 
eigenmode reverses its direction of propagation. We expect that 
friction-induced instabilities are also relevant in other astrophys- 
ical contexts, such as the angular momentum transport in accre- 
tion discs. 

Acknowledgements. The author thanks Manfred Schiissler, Robert Cameron, 
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Appendix A: Linear perturbations 

We summarise the expressions describing adiabatic perturba- 
tions of toroidal magnetic flux tubes in mechanical equilibriurn. 
Thei r detailed derivation can be found, for example, in ISchmitj 
([1998) and Fen-iz-Mas & Schussler (.1995,) . 

Small displacements, ^(s,t) = (^i,^„,^h), of the flux tube 
from its equilibrium position change the tangential, normal, and 
binormal vectors, t, n, b, respectively, of the Frenet basis, the 
curvature, k, the flow velocity, u, and the acceleration of the in- 
ternal plasma: 
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Quantities with index '0' and '1' are equilibrium and perturbed 
values, respectively, and indices \f and describe partial 
derivatives with respect to time, f, and equilibrium arclength, io- 
The flux tube is in pressure equilibrium with its environment. 
Assuming adiabatic perturbations, the field strength, density, and 
density contrast, Ap - p - p^, change according to the following 
equation: 
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with j6 = %npo/Bl and 
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The superadiabaticity 6 -V - Vad quantifies the stability of the 
stratification against convective turnover {5 > 0: convectively 
unstable). From the integrated Walen equation, ^ - ^ff^^ fol- 
lows 
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Appendix B: Coefficients of dispersion polynomials 

B.1. General case 

Based on the linearised equations of motion in the form of Eq. 
$Tl\ . the eigenfrequencies of magnetic flux rings in mechanical 
equilibrium parallel to the equatorial plane are determined by the 
monic dispersion polynomial 
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which stems from the coefficient matrix of Eq. ( |27] l. The coeffi- 
cients are: 

C5 — —6mRo + 2ia (B.2) 
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B.2. Axi-symmetric perturbations 

In the special case of axial symmetric perturbations the disper- 
sion equation reduces to the monic quartic polynomial in Eq. 
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( [30] l. The friction-independent part of the coefficients are: 
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C2 = 2 - 2 («o + 2) - 4 - ^1 w^gv (B.8) 



Co 



1 \^ 

4a;^BV (^o + 1)' - ( — I Ro" (Ro + if 



r/o 



B.3. Flux tubes in the equatorial piane 

For flux rings in the equatorial plane the dispersion equation fac- 
torises in a quadratic and quartic polynomial, which describe the 
perturbation perpendicular and parallel to the equatorial plane, 
respectively. The coefficients of the latter are: 

C3 = -AmRo + ia (B.IO) 
C2 = 4 [nf - l) Ro^ - 4 (mH 2) /Jo - 4 
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